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LENSING AND THE CENTERS OF DISTANT EARLY-TYPE GALAXIES 



(N 
O 
O 
(N 

C 
3 



> 
CO 

CN 

O 
CN 
O 

^ ' 

Or 
i 1 

O ■ 

a: 



Charles R. Keeton 1 

Astronomy and Astrophysics Department, University of Chicago, 5640 S. Ellis Ave., Chicago, IL 60637 

Submitted to The Astrophysical Journal 

ABSTRACT 

Gravitational lensing provides a unique probe of the inner 10-1000 pc of distant galaxies (z ~ 0.2-1). 
Lens theory predicts that every strong lens system should have a faint image near the center of the 
lens galaxy, which should be visible in radio lenses but have not been observed. We study these "core" 
images using models derived from the stellar distributions in nearby early- type galaxies. We find that 
realistic galaxies predict a remarkably wide range of core images, with lensing magnifications spanning 
some six orders of magnitude. More concentrated galaxies produce fainter core images, although not 
with any simple, quantitative, model independent relation. Some real galaxies have diffuse cores and 
predict bright core images (magnification ^/ C orc ^ Q.l) 5 but more common are galaxies that predict faint 
core images (/i CO ro 0.001). Thus, stellar mass distributions alone are probably concentrated enough 
to explain the lack of observed core images, and may require observational sensitivity to improve by 
an order of magnitude before detections of core images become common. Two-image lenses will tend 
to have brighter core images than four-image lenses, so they will be the better targets for finding core 
images and exploiting these tools for studying the central mass distributions of distant galaxies. 



Subject headings: galaxies: elliptical and lenticular, cD 
gravitational lensing 



galaxies: nuclei — galaxies: structure 



1. INTRODUCTION 

Galaxy centers are interesting places to study dynamics 
and galaxy formation. Their short crossing times make 
them sensitive to dynamical processes such as relaxation 
and binary black hole heating (e.g., Ebisuzaki, Makino 
& Okumura 1991; Milosavljevic & Merritt 2001). Their 
deep potential wells collect remnants of the galaxy forma- 
tion process such as the cores of accreted galaxies (e.g., de 
Zeeuw & Franx 1991; Barnes & Hernquist 1992). Their 
dark matter content provides clues to the interaction be- 
tween baryons and dark matter by adiabatic compression 
during galaxy formation (e.g., Blumenthal et al. 1986), 
and may even reveal properties of the dark matter parti- 
cle such as cross sections for self- interactions (e.g., Spergel 
& Steinhardt 2000). 

Nearby, galaxy centers can be studied directly with high 
spatial resolution observations. Hubble Space Telescope 
imaging of early- type galaxies shows that, contrary to the- 
oretical expectations (e.g., Tremaine 1997), the luminos- 
ity profiles diverge at small radii (e.g., Faber et al. 1997; 
Ravindranath et al. 2001; Rest et al. 2001). The profiles 
seem to fall into two classes: "core" galaxies have a dis- 
tinct transition between a steep outer profile and a shallow 
inner core that has I oc i?~ 7 with 7 < 0.3; while "power 
law" galaxies show no such break and have steep central 
cusps with 7 > 0.5. Interestingly, the global properties of 
the galaxies seem to correlate well with the centers. Core 
galaxies tend to be luminous, slowly rotating systems with 
boxy or elliptical isophotes, while power law galaxies tend 
to be faint, rapidly rotating systems with disky isophotes. 
Although the division may not be as sharp as originally 
thought (see Ravindranath et al. 2001; Rest et al. 2001), 
it still puts strong constraints on the formation process. In 
hierarchical merging scenarios, simple models cannot eas- 
ily explain why the large galaxies are so much less dense 
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than their small progenitors (see §7 of Faber et al. 1997, 
and references therein), and some additional process such 
as heating by binary black holes may be required (e.g., 
Milosavljevic & Merritt 2001; Milosavljevic et al. 2001). 

For distant galaxies, we cannot directly resolve 10-100 
parsec scales, but we can instead turn to a unique indi- 
rect probe offered by gravitational lensing. Lens theory 
predicts that if the central mass distribution is shallower 
than p oc r~ 2 then any multiply- imaged gravitational lens 
must have an odd number of images (Burke 1981; Schnei- 
der, Ehlers & Falco 1992). Standard image configurations 2 
have two or four bright images lying ~3-10 kpc from the 
center of the lens galaxy, with the remaining image just 10- 
100 pc from the center and much fainter than the others. 
Because the "core" image is very sensitive to the central 
surface density of the lens galaxy, with a higher density 
corresponding to a fainter image, it offers a unique way 
to constrain the density on scales that cannot be directly 
resolved. This probe of galaxy centers can in principle be 
applied to all lens galaxies, which now number more than 
60 and are predominantly early-type galaxies spanning the 
redshift range z ~ 0.3-1 (e.g., Kochanek et al. 2000). It 
is conceptually equivalent to using radial arcs to constrain 
the central profiles of lensing clusters (e.g., Mcllicr, Fort & 
Kneib 1993; Smail et al. 1996; Molikawa & Hattori 2001; 
Oguri, Taruya & Suto 2001). 

The best observational data on core images come from 
radio lenses, because the lack of radio emission from most 
lens galaxies enables sensitive searches for core images. 
The Cosmic Lens All-Sky Survey found 18 radio lenses 
but no core images, based on radio maps where the dy- 

2 The rare exception is a configuration where the source lies in 
a naked cusp (e.g., Schneider et al. 1992); among more than 60 
known lenses, APM 08279+5255 is the only candidate naked cusp 
lens (Lewis et al. 2002). The only other exception is B1359+154, a 
unique lens where three lens galaxies jointly produce six bright im- 
ages, and models predict three additional core images (Rusin et al. 
2001). 
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namic range is typically several tens to several hundreds 
but reaches 1200 for B1030+074 and 2000 for B0218+357 
(Rusin & Ma 2001; Norbury et al. 2002). Several other 
radio lenses have candidate core images (MG 1131+0456, 
Chen & Hewitt 1993; PMN J1632-0033, Winn et al. 2002) 
although the hypothesis that the central radio flux origi- 
nates in the lens galaxy cannot be ruled out. At optical 
and near-infrared wavelengths, APM 08279+5255 has an 
odd number of images (Ibata et al. 1999; Egami et al. 
2000), but its interpretation is not clear. The third image 
may be a core image, in which case it indicates a large low- 
density core in the lens galaxy (Ibata et al. 1999; Egami 
et al. 2000; Muhoz, Kochanek & Keeton 2001), or it may 
be a case of a "naked cusp" image configuration, in which 
case it contains no information about the center of the 
lens galaxy (Lewis et al. 2002). No other optical core im- 
ages have been seen, although the searches are of course 
hindered by light from the lens galaxies. 

The apparent discrepancy between data and theory pro- 
vides the desired opportunity to learn about the centers of 
distant galaxies. Motivated both by theoretical expecta- 
tions (see Tremaine 1997) and by ease of use, many anal- 
yses have assumed models with a finite density core and 
obtained limits on lens galaxy core radii (e.g., Narayan, 
Blandford & Natyananda 1984; Narasimha, Subramanian 
& Chitre 1986; Blandford & Kochanek 1987; Hinshaw & 
Krauss 1987; Narayan & Schneider 1990; Wallington & 
Narayan 1993; Kochanek 1996; Evans & Hunter 2002). 
Rusin & Ma (2001) instead used power law models and 
obtained a lower limit on the power law index, 7 > 0.8 
at 90% confidence for a surface density £ cx i?~ 7 . The 
question remained, though, whether these two classes of 
models were realistic enough to provide robust, model in- 
dependent conclusions about lens galaxy centers. Munoz 
et al. (2001) introduced double power law models where 
the core region is allowed to have a power law cusp whose 
index is independent of the density profile at large radii. 
They found that the lack of a core image in B1933+503 
robustly implies 7 > 0.6 for that one galaxy. Keeton 
(2001) studied core images statistically using models with 
cuspy stellar components, treated as generalized Hcrnquist 
(1990) models, embedded in dark matter halos. He found 
that the models were inconsistent with the data, perhaps 
because generalized Hernquist models may not accurately 
represent the stellar components of galaxies on 10-100 pc 
scales. 

The goal of this paper is to reconsider the core image 
problem using more realistic models derived from nearby 
galaxies, and more generally to discuss using core images 
as tools for studying the centers of distant (z ~ 0.3-1) 
galaxies. Nearby early-type galaxies have surface bright- 
nesses that can be modeled as a Nukcr law (Lauer et al. 
1996; Byun et al. 1996), and in §2 we discuss lensing by 
such galaxies. In §3 we consider in a general way what 
physical properties of lens galaxies determine core images 
properties, or conversely what we can learn about lens 
galaxies by studying core images. In §4 we study in de- 
tail the core images predicted by a sample of realistic lens 
galaxies. Finally, in §5 we offer a summary and conclu- 
sions. We assume the popular ACDM cosmology with 
matter density CIm — 0.3, cosmological constant = 0.7, 
and Hubble constant H n — 75 km s" 1 Mpc -1 . 



2. NUKER LAW LENSES 

The lensing properties of a galaxy with projected mass 
density £ are given by the lensing potential <j> that satis- 
fies the two-dimensional Poisson equation V 2 </> = 2S/S cr . 
Here S cr = (c 2 D os ) / ' (AttGD \D 0S ) is the critical surface 
density for lensing, where D Q \, D os , and As are angular 
diameter distances between the observer ("o"), the lens 
("1"), and the source ("s"). We consider a fiducial lens- 
ing situation with a lens galaxy at redshift z\ = 0.5 and a 
source at redshift z s = 2, which yields a critical density of 
S cr = 2230 Mq pc~ 2 for our adopted ACDM cosmology. 
The lensing deflection is given by a — V0, and the mag- 
nification depends on the second derivatives of <p. See the 
book by Schneider et al. (1992) for a full discussion of lens 
theory. 

Because lensing selects galaxies by mass, the sample of 
observed lens galaxies is dominated by early-type galaxies. 
In many nearby early-type galaxies the surface brightness 
distribution is well described by a Nuker law (Lauer et al. 
1995; Byun et al. 1996), 
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where 7 and (3 are the inner and outer power law indices, 
respectively, is the radius where the break in the power 
law occurs, a gives the sharpness of the break, and R is the 
surface brightness at the break radius. If the luminosity 
distribution has circular symmetry and the mass-to-light 
ratio is T, the lensing deflection is 
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where Kb = T/f,/E cr is the surface mass density at the 
break radius in units of the critical density for lensing, 
and 2*1 is the hypergeomctric function. If the stellar dis- 
tribution has ellipsoidal symmetry, the lensing deflection 
must be computed numerically using the formalism given 
by, e.g., Schramm (1990). 

Most galaxies contain central, supcrmassive black holes 
(e.g., Magorrian et al. 1998), so we consider adding them 
to the model. The deflection from a black hole is abh(R) = 
R\ I R where the black hole's Einstein radius is (in angular 
units) 
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We normalize the black holes using the observed corre- 
lation between the black hole mass M. and the velocity 
dispersion a of the parent galaxy (Gebhardt et al. 2000; 
Merritt & Ferrarese 2001). The net deflection is simply 
the sum of the deflections from the Nuker component and 
the black hole. 

As an example we consider the nearby galaxy NGC4486, 
which has Nuker parameters a = 2.82, (i = 1.39, and 
7 = 0.25 and is a fairly typical (albeit massive) exam- 
ple of the early-type galaxies in the sample studied by 
Faber et al. (1997). We imagine a lens galaxy obtained 
by moving NGC4486 to a typical lens redshift zi — 0.5. 
The deflection profile a(R) for the resulting lens is shown 
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Fig. 1. — The deflection profile for a mock lens galaxy obtained 
by taking the nearby galaxy NGC4486 and moving it to redshift 
Z[ = 0.5. The vertical dotted line indicates the Nuker break radius 
rj,. The dots mark the critical radii i? e in and R r ad- The top panel 
shows the result for the Nuker galaxy alone, while the bottom panel 
adds a central black hole normalized by the M.-cr relation from 
Gcbhardt et al. (2000). The mean core image magnification (^t C orc) 
is computed with the method presented in §3. 
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Fig. 2. — Cumulative distributions of magnification factors for the 
core images predicted by the NGC4486 mock lens galaxy, (a) Total 
distributions plotted for a circular (solid line) or flattened (dashed 
line) lens galaxy, (b) Distributions plotted separately for lenses 
with two or four bright images (doubles or quads), for a lens galaxy 
with ellipticity e = 0.5. There is more statistical uncertainty in the 
curve for quads than for doubles, because of the 7917 random source 
positions we examined only 238 of them corresponded to quads. 



in Fig. 1. The asymptotic behavior is a(R) oc i? 1 ~ 7 for 
R <C rb (if there is no black hole), and a(R) oc R}~@ for 
R TV Typical galaxies have 7 < 1 and (3 > 1, so the 
deflection is zero at the origin, rises to some finite peak, 
then slowly declines. There are two important radii cor- 
responding to the "critical curves," or curves along which 
the lensing magnification is infinite. The tangential crit- 
ical curve lies at the Einstein radius i? e in, which is the 
solution to a(i? j n ) = i? i n . The radial critical curve lies 
at the radius i? rac j which is the solution to da/dR = 1. 
The radial critical curve maps to a caustic at 

Umax = "(-Rrad) ~ -Rrad ■ (4) 

This caustic bounds the multiply-imaged region; sources 
with u < u m a X are multiply- imaged, while sources with 
u > w max are not. The core images are always contained 
in the region bounded by the radial critical curve. The 
radii i? e j n and i? rac j are marked in Fig. 1. 

If the galaxy has a steep central cusp with 7 > 1, the 
deflection diverges at the origin and is a monotonically 
decreasing function of radius. In this case, the radial crit- 
ical curve does not exist and the lens never produces a 
core image. If the galaxy contains a central black hole, 
the black hole causes the deflection to diverge at the ori- 
gin and suppresses some of the core images (Mao, Witt 
& Koopmans 2001). However, the black hole changes the 
deflection only at very small radii, so it has little effect on 
the critical radii i? ; n and i? r ad or on the mean core image 
magnification (/i C orc) (see Fig. 1). Adding a dark matter 
halo to the lens galaxy would raise the outer deflection 
profile and make it approximately flat, but it would have 
little effect on the central deflection profile that determines 
the properties of the core images unless it were much more 



centrally concentrated than the light. 

To characterize the core images expected in this lens, we 
study the core image magnification distribution. The dis- 
tribution can be computed fairly rapidly using inverse ray 
shooting (e.g., Kayser, Refsdal & Stabell 1986; Wambs- 
ganss 1997). Fig. 2a shows the distribution for circular and 
flattened lens galaxies. The distribution is broad, span- 
ning more than two orders of magnitude, with a median 
value /i coro = 0.0019 and a mean value (fJ, COTe ) = 0.0053 
for the circular case. Making the galaxy flattened shifts 
the magnification distribution to higher values; but even 
for ellipticity e = 0.5 the shift is only 0.05 dex, so the core 
image magnification distribution is largely insensitive to 
ellipticity in the lens galaxy. (We would see more of an ef- 
fect if we described core images by their flux ratio relative 
to the bright images, /i C oro/Mbright, because /^bright is quite 
sensitive to ellipticity.) 

Ellipticity is important in one respect. If the lens galaxy 
is non-spherical, some source positions correspond to lenses 
with two bright images (doubles) , while others correspond 
to lenses with four bright images (quads). (See, e.g., Schnei- 
der et al. 1992.) Fig. 2b shows the core image magnifica- 
tion distributions for quads and doubles separately. Be- 
cause the source lie closer to the origin for quads than for 
doubles, and the core image magnification increases with 
the distance of the source from the origin, quads tend to 
have smaller core image magnifications than doubles. The 
quantitative details depend on the ellipticity, on any exter- 
nal tidal perturbation (shear) that may affect the lens, and 
on the stellar profile of the galaxy, but the general result 
is robust: quads to tend have smaller /i coro than doubles. 

The ellipticity effect combines with an observational se- 
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lection effect. Quads generally have larger total magnifi- 
cations than doubles, so the sources in quads tend to be 
intrinsically fainter than the sources in doubles. Together, 
their fainter sources and smaller core image magnifications 
suggest that quads will tend to have fainter core images 
than doubles. The lack of core images in observed quads 
should therefore be less surprising than in doubles. Con- 
versely, doubles should be better systems than quads for 
searching for core images. 

Thus, while ellipticity in the lens galaxy is important 
in understanding differences between quads and doubles, 
it is not very important in the overall distribution of core 
image magnifications. In the remainder of the paper we 
therefore neglect ellipticity and use circular lens galaxies. 

3. WHAT DO CORE IMAGES PROBE? 

In order to derive meaningful conclusions from observa- 
tional constraints on core images, it is important to under- 
stand how core images depend on the physical properties 
of lens galaxies. (Merely understanding the parameter de- 
pendencies in parametric lens models does not allow strong 
physical conclusions.) Although the full distributions of 
core images predicted for a particular lens galaxy must be 
computed numerically, the mean magnification (p C orc) can 
be studied analytically. Moreover, in §4 we argue that this 
quantity is actually a good way to characterize the distri- 
bution. Hence, in this section we study (p C orc) in general 
terms. 

Formally, we have 

/mult Mcore(u) du 
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f ,, du 

J mult 



f dx 

J core 
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The first line is simply the definition of the average, where 
the integral extends over the multiply- imaged region of 
the source plane. In the second line, the integral in the 
numerator extends over the "core" region in the image 
plane, defined to be the region within the radial critical 
curve; this equality holds because ^ CO rc = \dx/du\ is the 
Jacobian of the transformation between the source and 
image planes. In words, eq. (6) says that the mean core 
image magnification is equal to the area within the radial 
critical curve divided by the lensing cross section (the area 
of the multiply-imaged region in the source plane). For 
circularly symmetric lenses, 

(/^corc) — (-^rad/^max) • (7) 

For non-circular lenses, Fig. 2 suggests that this is still a 
good approximation, because the ^ CO rc distribution is not 
terribly sensitive to ellipticity. These relations were first 
given by Keeton (2001). 

Now we focus on circularly symmetric lenses. Given the 
definition of u max from eq. (4), we can write eq. (7) as 

(^corc) = K-RradV-Rrad - 1]~ 2 , (8) 

= «^ ld -ir 2 , (9) 

where {h)r is the mean surface density within radius R, 
in units of the critical density for lensing. The second 
equality holds because by the definition of the deflection, 

a{R) * rR 
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Alternatively, returning to eq. (8) and using the identities 

2k(R), (11) 

= 1, (12) 
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we obtain 
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Eqs. (9) and (13) indicate that the mean core image mag- 
nification is given very simply from either the surface mass 
density at the radial critical curve i? ra( j or the mean surface 
mass density within i? rac j. 

What remains is to understand what physical proper- 
ties of the galaxy determine the radial critical curve i? ra d- 
Equating eqs. (9) and (13), we find that i? ra d is the solu- 
tion to 

2 K ( J R rad ) = («) iJrad + l. (14) 

This relation implies that the radial critical curve is related 
to the concentration of galaxy's (projected) mass distri- 
bution. We expect k(R) to be a decreasing function, so 
{k)r > k(R) for all R. If n{R) is steep (the mass is con- 
centrated), then (k) is large and we must go to large k 
(small R) to satisfy eq. (14); eq. (13) then implies faint 
core images. 3 Conversely, if k(R) is shallow then eq. (14) 
is satisfied at smaller k (larger R) , and the core images are 
brighter. 

These results can be demonstrated with two simple ex- 
amples. First, consider a softened isothermal sphere with 
surface mass density k(R) = (b/2)(s 2 + R 2 )- 1 / 2 , where 
s is a core radius and & is a scale radius that equals the 
Einstein radius in the case s = 0. The critical radii and 
mean core image magnification are 

R cin = [b(b - 2s)} 1/2 , 
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where £ = (46 + s) 1 / 2 and ( = s 1 / 2 . Decreasing the core 
radius reduces (p C orc), with (^ coro ) « s/b for s<6. Sec- 
ond, consider a power law density p oc r~ 7 or k oc i? 1-7 , 
with 7 > 1 to ensure that k(R) is a decreasing function. 
The critical radii and mean core image magnification are 

i(2 _ 7 )i/(7-i) ; (19) 



Rr&d 
^max 

(/^corc) 
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Increasing 7 (making the profile steeper) decreases (/i CO re)- 
These expressions are valid only for 1 < 7 < 2, because 
for 7 > 2 the density profile is so steep that the radial 

3 This result also explains how adding a central black hole affects 
the mean core image magnification. The black hole suppresses some 
core images, reducing (^t C oro) (Mao ct al. 2001). In the language of 
our analysis, the black hole increases without changing k(R), 

so we must move to larger k (smaller R) to keep eq. (14) satisfied. 



5 



critical curve does not exist and the model does not pro- 
duce core images. Previous studies used models like these 
to understand the inverse relation between the brightness 
of core images and the concentration of the lens galaxy, 
but eqs. (13) and (14) now give it in a general, model 
independent form. 

4. CORE IMAGES IN REALISTIC GALAXIES 

4.1. The sample 

To understand core images in realistic galaxies, we study 
models constructed from a sample of observed galaxies. 
This approach ensures that we examine the region of pa- 
rameter space occupied by real systems. We seek galaxies 
with well resolved luminosity profiles, plus measured ve- 
locity dispersions as mass indicators. Faber et al. (1997), 
Carollo et al. (1997), Carollo & Stiavelli (1998), and Ravin- 
dranath (2001) have published samples of nearby early- 
type galaxies observed with high-resolution Hubble Space 
Telescope imaging at optical or near-infrared wavelengths. 4 
Together the samples comprise 73 distinct galaxies with 
both Nuker law fits and velocity dispersions, which are 
summarized in Table 1. The different samples use differ- 
ent passbands: V for the Faber sample, H for the Ravin- 
dranath sample, and V and I for the Carollo sample. How- 
ever, we can check for wavelength dependence and other 
systematic effects because some of the galaxies appear in 
more than one sample: six galaxies in both the Faber and 
Ravindranath samples, five galaxies in both the Faber and 
Carollo samples, six galaxies in both the Ravindranath and 
Carollo samples, and ten galaxies with multiple passbands 
in the Carollo sample. (No galaxies appear in all three 
samples.) 

We use this sample to construct mock lens galaxies by 
moving each galaxy to a typical lens redshift z; = 0.5, 
and assuming a typical source redshift z s = 2. Following 
Faber et al. (1997), we compute the mass-to-light ratio 
for each galaxy using a spherical and isotropic dynami- 
cal model. Faber et al. (1997) give mass-to-light ratios 
for their sample, providing a check for our calculations. 
We then compute the lensing properties of the mock lens 
galaxies, which are summarized in Table 1. 

As discussed in §2 we focus on circularly symmetric 
galaxies because ellipticity has little effect on the over- 
all distribution of core image magnifications. We consider 
only the Nuker components of the galaxies, neglecting any 
nuclear components that may be indicated by fits to the 
surface brightness distribution. This is a conservative ap- 
proach to our problem, because any additional mass con- 
centration would decrease the predicted core image mag- 
nifications and bring the models closer to agreement with 
the data. We also neglect dark matter halos. Dark matter 
does not appear to be dynamically important in the inner 
5-10 kpc of elliptical galaxies (e.g., Gerhard et al. 2001). 
While it is important for lensing (because lensing depends 
on the projected mass; e.g., Treu & Koopmans 2002), dark 
matter would have little effect on the projected mass distri- 
butions on the <200 pc scales important for core images 
unless it were substantially more concentrated than the 
light. 

We can check that our mock lenses are reasonable in sev- 

4 Rest et al. (2001) give a similar sample, but without velocity dis- 
persions. 
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Fig. 3. — Relation between the velocity dispersion a and Einstein 
radius R c i n for the mock lens galaxies. Crosses and boxes indicate 
galaxies in the Faber and Ravindranath samples, respectively. The 
dotted line shows the relation i? e i n oc a 2 ; the fitted zero point is 
—4.84 (in log units), compared with a predicted value of —4.74 for 
comparable Singular Isothermal Spheres. 



eral ways. First, we compute the Einstein radius for each 
galaxy and compare it to the velocity dispersion in Fig. 3; 
this tests whether the lensing and dynamical masses are 
consistent. For comparison, a simple Singular Isothermal 
Sphere (SIS) lens has the scaling i? i n cx a 2 , with a zero 
point of —4.74 (in log units) for zi = 0.5 and z s = 2 in 
our adopted ACDM cosmology (see Schneider et al. 1992). 
The mock galaxies are consistent with this scaling and a 
fitted zero point of —4.84, although with a scatter of 0.16 
dex. The fact that the mock galaxies lie on average 0.1 
dex below the expected SIS relation may be related to 
our neglect of dark matter halos. Still, the dynamical and 
lensing properties are related in a sensible way suggesting 
that the mock lens galaxies are not unreasonable. 

As a second check, we consider the galaxies that ap- 
pear in more than one of the original samples. For ex- 
ample, for each galaxy that was observed and modeled by 
both Faber et al. (1997) and Ravindranath et al. (2001), 
we compute the lensing properties using both models and 
compare them. In this way we test whether the use of dif- 
ferent modeling techniques and different passbands affects 
our lensing results. Fig. 4 compares the Einstein radii and 
mean core image magnifications for all of the duplicate 
galaxies. There is fair agreement in the Einstein radii, 
although with some scatter because different studies find 
somewhat different values of the outer slope /?; the main 
outlier is NGC524, where dust is known to affect the lu- 
minosity profile at optical wavelengths (Lauer et al. 1995; 
Ravindranath et al. 2001). There is good agreement in the 
mean core image magnification, so our conclusions about 
the properties of core images are not sensitive to whose 
data we use. For the remainder of the paper, we adopt as 
our main sample all of the Ravindranath galaxies plus the 
Faber galaxies that are not in the Ravindranath sample. 
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4.2. A plethora of core images 

Fig. 5a shows the core image magnification distributions 
for nine of the mock lens galaxies to illustrate the range 
of effects. Each mock lens galaxy has a distribution of 
core image magnifications that spans some two orders of 
magnitude, and the complete set of mock lenses spans six 
orders of magnitude in ^ corc . Realistic lens galaxies predict 
a remarkably wide range of core images. 

Studying the full /x cor e distribution for each mock lens 
galaxy is impractical, so we would like to describe each 
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Fig. 5. — Cumulative core image magnification distributions for 
nine of the mock lens galaxies: the three with the lowest non-zero 
values for ((i corc ) (dotted lines; NGC221, NGC3377, NGC4621); the 
three with values for (//core) at the median of the sample (solid lines; 
NGC4570, NGC5982, NGC4649); and the three with the largest val- 
ues for (//core) (dashed lines; NGC4239, NGC6166, NGC5273). (a) 
Distributions of the core image magnification ficoic- (b) Distribu- 
tions of the normalized magnification /i cor c / (/icorc ) ■ 



galaxy by a single characteristic quantity. We propose 
to use the mean core image magnification (^ CO rc) as a 
good characteristic value, partly because this quantity is 
easy to compute (see §3), and partly because Fig. 5b sug- 
gests that it is indeed characteristic of the overall distribu- 
tion. Specifically, when plotted in terms of the normalized 
quantity ^ coro / (^ C oro) the distributions for a wide range of 
galaxies all lie on top of each other. Although there are 
differences in the shapes of the distributions at the faint 
end, the distributions at the bright end are remarkably 
similar. Hence, we believe that (/i CO ro) is a useful charac- 
terization of the distribution of core image magnifications 
for a particular galaxy, especially at the bright end. 

Fig. 6 shows a histogram of the (^ CO rc) values for the 73 
mock lens galaxies. (The values are given in Table 1.) 
Note that the histogram is intended only to show the 
broad range of core image properties in our sample; it 
should not be interpreted as a global distribution of (^ C orc) 
values, because our sample is not a statistical sample of 
galaxies. Still, the histogram is instructive in illustrating 
that the mean values (/x cor c) span more than four orders 
of magnitude, from (^ C orc) = 1-8 x 10~ 5 for NGC221 to 
(A*corc) = 0.19 for NGC5273 — plus two galaxies that never 
produce core images (NGC1172 and NGC4742). Again we 
see the wide range of core image properties in realistic lens 
galaxies; some lenses should have bright, detectable core 
images, while others should have core images that are es- 
sentially invisible. 

In principle the central supermassive black holes that are 
common in galaxies can suppress core images (Mao et al. 
2001), but in practice they have little effect. Table 1 gives 
the values for (^ C orc) when the galaxies contain black holes 
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Fig. 6. — Histogram of the mean core image magnifications for 
the 73 mock lens galaxies. Galaxies that do not produce core images 
are arbitrarily placed at log(/^ coro ) = — 6. 
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holes, for four of the galaxies from Fig. 5. Black holes can suppress 
core images (see Mao ct al. 2001), but only at the faint end of the 
core image distribution. 



normalized by the M.-cr relations measured by Gebhardt 
et al. (2000) and Merritt & Ferrarese (2001). In a few 
cases the black hole modifies the central potential enough 
to erase all core images (NGC221, NGC3115, NGC3377, 
NGC3900, NGC4464, NGC4467, NGC4621, and NGC5838). 
But all of these cases have (^ coro ) < 0.001 even without 
a black hole, which suggests that adding a black hole can 
erase all core images only if the core images are faint to 
begin with. In the remaining cases the black hole reduces 
(Mcore) by only about 0.1 dex. The black hole does sup- 
press some core images, but only at the faint end of the 
distribution, as shown in Fig. 7. In other words, black 
holes have little effect on the core images in realistic lens 
galaxies, especially at the bright end of the core image 
distribution. Black holes therefore fail to resolve the core 
image paradox. These conclusions are based on local mea- 
surements of the M,-a relation, which may not hold at 
higher redshifts; but black holes would have to be sub- 
stantially more massive (relative to their parent galaxies) 
at z ~ 0.5 than at z — in order to change our conclu- 
sions. Another possibility is that binary black holes arc 




Fig. 8. — Histograms of the lensing critical radius -R ra d, which 
bounds the core image region, (a) The critical radius in parsecs. 
(b) The ratio of the critical radius to the Nuker break radius rj,. 
Galaxies that lack a critical radius i? r ad are not included. 



common in galaxies at z ~ 0.5. In principle binary black 
holes can be more effective than single black holes at sup- 
pressing core images, but for realistic binary parameters 
the differences are small and binary black holes do not 
strongly suppress core images (Keeton & Zhao, in prep.). 
In the remainder of the paper we neglect black holes. 

We now seek to understand the diversity of core image 
properties in the mock lens galaxies, and to identify the 
galaxy properties that affect the core images. First, we 
examine the region of the galaxy that is probed by core 
images. Recall that the core images are always contained 
within the inner lensing critical curve, which has radius 
-Rrad (see §2). Fig. 8 shows that the core image region has 
a characteristic scale ~200 pc (the median value of i? ra( j 
for the sample), and is comparable to or slightly larger 
than the Nuker break radius. Thus, core images probe the 
inner tens to hundreds of parsecs at the centers of galaxies. 

Next, we consider how the mean core image magnifica- 
tion (^ coro ) depends on various properties of the galaxies, 
as shown in Fig. 9. The interesting general result is that no 
single property of a galaxy strongly determines its core im- 
age properties. There is no simple relation between (/z CO re) 
and the outer and inner Nuker power law indices (3 and 7 
(Figs. 9a and 9b). This result makes sense in combina- 
tion with Fig. 8: the lensing critical radius i? ra d is often 
comparable to the Nuker break radius rj,, so we are in a 
regime where the galaxy cannot be described as a simple 
power law, and (^ C orc) is not dominated by either (3 or 
7 separately. The distinction between "core" and "power 
law" galaxies seen in their luminosity profiles and dynam- 
ics (Faber ct al. 1997) does not appear to carry over into 
lensing and core images. There is likewise no simple rela- 
tion between (^ C oro) and galaxy mass (as represented by 
velocity dispersion; Fig. 9c); galaxies with a given a have 
(Mcorc) values that range over some three orders of magni- 
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Fig. 9. — Scatter plots of the mean core image magnifications (/i C orc) versus various galaxy properties. Crosses and boxes indicate galaxies 
in the Faber and Ravindranath samples, respectively. Galaxies that do not produce core images are arbitrarily placed at log(/^ C orc) = —6. 



tudc or more. Dwarf galaxies are not systematically more 
or less likely than giant galaxies to produce bright core 
images. 

Based on the arguments in §3 we expect a connection 
between the core image properties and some measure of 
the concentration of the mass distribution. There is a 
trend between (/i CO rc) and the surface density at the break 
radius (k^, Fig. 9d): galaxies with higher break densities 
tend to predict fainter core images. There is also a general 
trend between with the Nuker break radius (Fig. 9d): 
galaxies with smaller break radii, and hence smaller core 
regions, tend to predict fainter core images. But this trend 
is neither strong nor tight, because by itself cannot dis- 
tinguish between galaxies that are large and highly con- 
centrated and those that are less concentrated but intrin- 
sically small. One way to remove this effect is to normalize 
the break radius using the Einstein radius as a measure of 
the global scale of the galaxy; because i? i n oc cr 2 , we ac- 
tually use rfc/cr 2 (Fig. 9f). These three trends all indicate 
that there is indeed some connection between the core im- 
ages and the concentration of the galaxy such that more 
concentrated galaxies tend to predict fainter core images. 
However, the trends have significant scatter and thus are 
not highly predictive. In other words, there does not ap- 
pear to be a simple measure of a galaxy's concentration 
that strongly determines its core image properties. 



To summarize, realistic lens galaxies have an extremely 
wide range of core image properties; some should have 
bright, detectable core images, while others should have 
core images that are very faint or absent altogether. A 
galaxy's core image properties are related to its mass con- 
centration, but there does not appear to be any simple 
measure of concentration that yields a clean prediction of 
the core image properties. The complication is that in 
Nuker law lenses the lensing critical radius i? ra( j tends to 
be comparable to the break radius r^, which means that 
all of the Nuker parameters affect the core image proper- 
ties. This conclusion has a somewhat surprising corollary. 
We might expect that the distinction between "core" and 
"power law" galaxies, with their different luminosity pro- 
files, would be obvious in their core image properties, but 
it is not. 

4.3. Should we see core images? 

We can now re-evaluate the core image problem in terms 
of our expectations for realistic lens galaxies. Although our 
sample is not a proper statistical sample of galaxies, if we 
assume that it does at least represent the range of realistic 
galaxy properties then we can consider whether the lack 
of observed core images is surprising or not. The main 
issue is to understand what types of lens galaxies produce 
bright core images and whether we should expect to find 
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many lenses from such galaxies. There are three galaxies 
in the sample with (^ C orc) > 0.1: two are dwarf galax- 
ies (NGC4239 and NGC5273), while the third is a giant 
elliptical at the center of a cluster (NGC6166). Because 
lensing selects galaxies by mass it tends to select against 
dwarf galaxies, so even if some dwarf galaxies are good at 
producing core images they are unlikely to produce lenses 
in the first place. As for NGC6166, it is unusual for be- 
ing a brightest cluster galaxy as well as the most distant 
(120 Mpc) and least concentrated (lowest surface bright- 
ness and largest break radius) galaxy in the Faber et al. 
(1997) sample. Thus, it appears to be an atypical galaxy 
drawn from the tail of the galaxy population. Statisti- 
cally, lensing is not likely to select rare galaxies to be lens 
galaxies. 

Lensing selection effects are important for core images 
in another way. Consider a set of galaxies that are sim- 
ilar to each other, specifically a set of galaxies with the 
same Einstein radius. Within this set, are galaxies that 
predict bright core images any more or less likely to be 
selected for lensing than galaxies that predict faint core 
images? Fig. 10 shows the lensing cross section, normal- 
ized by the area within the Einstein ring, versus (/U C ore)- 
There is a clear decrease in the normalized cross section 
with increasing (^ CO re), and it does not depend on whether 
magnification bias is included or omitted. The trend for 
the mock lens galaxies agrees well with analytic predic- 
tions for simple softened power law lens models. At fixed 
Einstein radius, then, galaxies that predict bright core im- 
ages have smaller lensing cross sections than galaxies that 
predict faint core images, so they are less likely to be se- 
lected for lensing. There is an intrinsic bias against lenses 
with bright core images. 

Thus, the types of galaxies that can produce bright core 
images are probably not common, and they have small 
lensing cross sections relative to galaxies that are com- 
parably massive but produce faint core images. The two 
effects combine to suggest that bright core images are not 
likely to be prevalent in observed lens samples. 

4.4. Comparison with data 

The observational constraints on core images take the 
form of upper limits on the core image flux / cor c in ob- 
served lenses. The magnification cannot be directly con- 
strained because the intrinsic flux of the source is un- 
known. One way around this problem is to use the flux 
ratio of the core image to the brightest image, because the 
source flux factors out to leave f COIC / /bright = Mcorc/Mbright- 
This approach requires explicitly solving the lens equation 
to find all images, and thus is valuable for applications like 
lens modeling where the lens equation must be solved any- 
way (e.g., Muhoz et al. 2001). While it can be used for sta- 
tistical analyses (e.g., Rusin & Ma 2001; Keeton 2001), it 
eliminates technical simplifications like the ability to com- 
pute (^coro) rapidly without solving the lens equation (but 
see Evans & Hunter 2002 for a different kind of technical 
simplification). It also makes the conclusions sensitive to 
quantities like cllipticity in the lens galaxy, because /i CO rc 
is fairly insensitive to ellipticity but ^bright is not. 

An alternate approach is to fit a lens model to the ob- 
served images to constrain the source flux, and then com- 
bine the inferred source flux with the limits on / cor c to 
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Fig. 10. — Scatter plots of the lensing cross section versus the 
core image magnification. The cross sections are normalized by the 
area within the Einstein ring, ^R 2 in - (Top) A\ is the simple lensing 
cross section, or the area of the multiply-imaged region in the source 
plane. The curves show results for simple softened power law lens 
models S oc (s 2 +R 2 )-l 3 with /3 = 1.0, 1.1, 1.2, 1.3, 1.4, and 1.5 from 
bottom to top. (Bottom) A2 is the lensing cross section corrected 
for magnification bias. 



put limits on /i cor c- This approach should be robust be- 
cause the model properties that determine the source flux 
from the observed images (the enclosed mass on 3-10 kpc 
scales) decouple from the properties that affect core im- 
ages (the density profile on <200 pc scales). With this 
approach, Norbury et al. (2002) obtain upper limits on 
Hoove for 15 radio lenses from the Cosmic Lens All-Sky 
Survey, as shown in Fig. 11. 

For an accurate comparison we need a prediction of the 
overall /i C orc distribution from the models (not merely the 
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Fig. 11. — Upper limits on core image magnifications for 15 ob- 
served radio lenses (Norbury et al., in prep.). The points show 5a 
upper limits on ^ C ore- (a) The data for all lenses. The dashed curve 
shows the net /i C orc distribution predicted by a weighted sum of 
the model galaxies (see text), (b) The data for doubles and quads 
shown separately. 

set of (/Ucore) values for the sample galaxies, as in Fig. 6). 
A proper prediction is impossible because, again, our sam- 
ple is not a statistical sample of galaxies. Nevertheless, for 
the sake of comparison we naively combine our sample by 
summing all of the ^ C oro distributions (a few of which are 
shown in Fig. 5), weighting each galaxy by its lensing cross 
section. The result is shown as the dashed curve in Fig. 11. 
It is interesting to see that the model predictions lie well 
below the limits from the data. The model does not pre- 
dict that core images should be brighter than observed. 
In fact, it suggests that the observational sensitivity may 
need to improve by more than an order of magnitude be- 
fore detections of core images become common. Despite 
concerns that the model sample is not statistically com- 
plete, we believe that the general conclusion is reliable. 
Galaxies that produce bright core images would have to 
be substantially more common in the universe than in our 
sample in order to make the model predictions inconsistent 
with the current observational data. 

Evans & Hunter (2002) reach similar conclusions from 
an analysis of softened power law potential models. They 
argue that the break radii in the galaxies observed by 
Faber et al. (1997) are small enough to make core images 
faint even if the galaxies have simple finite-density cores. 
Although we believe our models to be more realistic be- 
cause they allow more general cores and are constructed 
directly from the fitted profiles, it is reassuring that the 
conclusions are consistent. 

We previously found that Hernquist model galaxies pre- 
dict core images to be more common than observed (Kee- 
ton 2001). Like a Nukcr law, a Hernquist model has a steep 
outer profile that smoothly changes to a shallow central 
cusp, but with a larger transition radius. Our new mod- 
els should be more realistic because they are based on fits 



to the luminosity profiles of observed galaxies, which have 
small transition radii. The smaller break radii mean higher 
central densities and thus fainter predicted core images. 

Another interesting result from the data is that the up- 
per limits from quads are weaker than the upper limits 
from doubles (Fig. lib). The flux limits for the various 
lenses are comparable, but the quads tend to have fainter 
sources. Combining this result with our prediction that 
quads will tend to have fainter core images than doubles 
(Fig. 2b), we conclude that the lack of core images in ob- 
served quad lenses is no surprise at all. 

5. CONCLUSIONS 

Core images in strong gravitational lens systems provide 
a unique probe of the centers of galaxies at rcdshifts z ~ 
0.2-1. The brightnesses of core images are determined by 
the density profiles of galaxies inside <200 pc. The lack of 
core images in observed lenses, especially in radio lenses, 
sets strong lower limits on the central densities of the lens 
galaxies (e.g., Wallington & Narayan 1993; Munoz ct al. 
2001). 

The mapping between core images and galaxy centers 
can be studied in two directions. In the forward prob- 
lem knowledge of galaxy centers is used to make predic- 
tions about core images. Based on the first lens models 
drawn directly from the resolved stellar mass distributions 
of nearby early-type galaxies, we predict that real galax- 
ies should produce a remarkably wide range of core im- 
ages. Some should have bright core images (magnification 
Mcore 0.1), while many others will have core images that 
are faint (^ cor o ^ 0.001) or absent altogether. 

Qualitatively, more concentrated galaxies produce fainter 
core images. Quantitatively, however, there does not seem 
to be a simple predictive relation between observed galaxy 
properties and core images. Lensing is biased against 
galaxies with bright core images, because they have smaller 
cross sections than comparable galaxies with faint core im- 
ages. Four-image lenses should tend to have fainter core 
images than two-image lenses, because in quads the source 
is always close to the center of the lens galaxy where the 
core image magnification is low. Supermassive black holes 
in the centers of galaxies can suppress faint core images, 
but they have little effect on bright core images or on the 
mean magnification. 

The connection between core images and galaxy cen- 
ters can also be studied in the inverse problem, where the 
analysis of core image data yields constraints on the cen- 
ters of lens galaxies. Previous work placed limits on the 
core radius or on the logarithmic slope of the density (e.g., 
Wallington & Narayan 1993; Rusin & Ma 2001), but it was 
not clear how model dependent those constraints were. We 
obtain a general statement of the connection between the 
density profile and core images: the mean core image mag- 
nification (^corc) is inversely related to the density at the 
lensing critical radius i? rac i (eq. 13); and this critical radius 
is determined by the shape of the density profile (eq. 14), 
with more concentrated galaxies corresponding to smaller 
critical radii and fainter core images. Unfortunately, nei- 
ther our general formalism nor our Nuker law lenses sug- 
gest any simple, model independent measure of the mass 
concentration that determines the core image properties. 
The interpretation of core image data will therefore con- 



tinue to rely on detailed models of individual lenses. The 
model dependence can be held in check, though, by us- 
ing general models like the Nuker law or the cuspy lenses 
introduced by Muhoz et al. (2001), as opposed to overly 
simple fiat core or pure power law models. 

We conclude that in many cases the stellar mass in lens 
galaxies is probably concentrated enough to render core 
images faint (also see Evans & Hunter 2002). This is not 
to say that bright core images cannot exist — certainly 
there are realistic galaxies that predict bright core images, 
and the probability that they are selected for lensing is 
non-zero. But the fact that core images have not yet been 
found (with perhaps one or two exceptions) is probably 
not a surprise. As the search continues, two-image lenses 
should be better targets than four-image lenses for reveal- 
ing core images. 
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Table 1 
Galaxy Sample 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Name 


D 


a 


a 





7 


log K b 


log r b 


log -Rein 


log -Rrad 


log(Mcoro) 


log(Alcoro) 


log(Alcoro) 


Ref 




(Mpc) 


(km/s) 










(arcscc) 


(arcscc) 


(arcsec) 











NGC221 


0.9 


85 





.98 


1 


.36 





,01 


2 


,86 


-4.02 


-1 


.26 


-3 


.17 


-4 


.67 


-6 


.00 


-6.00 


1 




0.7 


77 


4 


.66 


1 


.26 





.50 


2 


.56 


-3.71 


— 1 


.29 


—3 


.28 


-4 


.74 


-6 


.00 


-5.31 


2 


NGC224 


0.9 


220 


4 


.72 





.87 





12 


3 


,25 


-3.65 





.43 


— 


.58 


— 1 


.65 


— 1 


.66 


-1.66 


1 


NGC474 


32.5 


169 


1 


.23 


1 


.90 





37 





.91 


— 1.39 


— 


,43 


— 1 


.30 


— 2 


.08 


—2 


.10 


— 2.10 


2 


NGC524 


24.6 


275 


1 


.29 


1 


.00 





,00 


1 
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1 


.05 


1 


.43 





11 


1 


,94 


— 2.22 


-0 


26 


— 1 


.58 


—3 


.15 


—3 


.28 


-3.32 


2 


NGC1172 


31.8 


113 


1 


.52 


1 


.64 


1 


,01 





.18 


-1.21 


-0 


,76 


-6 


.00 


-6 


.00 


-6 


.00 


-6.00 


1 


NGC1316 


19.1 


240 


1 


.16 


1 


.00 





,00 


1 


.85 


— 2.21 





,20 


— 


.87 


— 1 


.96 


— 1 


.99 


— 2.01 


1 


NGC1399 


19.1 


333 


1 


.50 


1 


.68 





,07 


1 


.49 


-1.33 





10 


— 


.95 


—2 


.66 


-2 


.69 


— 2.72 


1 


NGC1400 


22.9 


265 


1 


.39 


1 


.32 





00 


2 


.08 


— 2.22 


-0 


09 


— 1 


49 


-3 


.29 


—3 


44 


-3.52 


1 


NGC1426 


22.9 


150 


3 


.62 


1 


.35 





.85 





.89 


-1.53 


-0 


.51 


— 1 


.52 


-2 


.25 


—2 


.31 


-2.30 


1 


NGC1600 


53.5 


340 


1 


.98 


1 


.50 





,08 


1 


.02 


-0.88 





26 


— 


.52 


-1 


.68 


— 1 


.69 


-1.70 


1 


NGC1700 


37.9 


230 





.90 


1 


.30 





,00 


2 


.23 


-2.57 


-0 


.19 


— 1 


.62 


-3 


.27 


-3 


.51 


-3.59 


1 




54.1 


230 





.46 


1 


.65 





,01 


1 


.80 


-2.08 


-0 


,17 


— 1 


.40 


—2 


.80 


—2 


.93 


-2.96 


3V 




54.1 


230 





.47 


1 


.68 





,01 


1 


.81 


-2.08 


-0 


19 


— 1 


.44 


—2 


.89 


-3 


.03 


-3.07 


31 


NGC2636 


35.7 


85 


1 


.84 


1 


.14 





04 


1 


,41 


-2.59 


-0 


,88 


— 1 


.87 


-1 


.97 


—2 


.01 


-1.99 


1 


NGC2685 


16.2 


114 


1 


.69 


1 


,52 





.73 





,64 


— 1.48 


-0 


68 


— 1 


.56 


-1 


87 


— 1 


.90 


-1.89 


2 


NGC2832 


96.2 


330 


1 


.84 


1 


40 





02 


1 


,27 


-1.16 





23 


-0 


.68 


-2 


.02 


—2 


.04 


-2.06 


1 


NGC2841 


14.1 


229 





.93 


1 


.02 





.01 


2 


.35 


— 2.84 





.08 


-1 


.17 


-2 


.41 


—2 


.53 


-2.56 


1 


NGC3115 


9.0 


280 


1 


.47 


1 


.43 





.78 


1 


.60 


-1.69 


-0 


.13 


-1 


.61 


—3 


.59 


-4 


.18 


-6.00 


1 




6.7 


264 


1 


.13 


1 


.80 





,73 


1 


.36 


-1.51 


-0 


,25 


— 1 


.65 


-3 


.62 


-4 


.00 


-6.00 


2 


NGC3377 


10.6 


152 


1 


.92 


1 


.33 





29 


2 


.53 


-3.12 


-0 


74 


—2 


.57 


-4 


.48 


-6 


.00 


-6.00 


1 


NGC3379 


10.6 


225 


1 


.59 


1 


.43 





,18 


1 


.61 


-1.84 


-0 


21 


— 1 


.35 


-2 


.72 


—2 


.78 


-2.80 


1 




8.1 


209 


1 


.82 


1 


.45 





18 


1 


.71 


-1.96 


-0 


29 


— 1 


.50 


—2 


.95 


—3 


02 


-3.03 


2 


NGC3384 


8.1 


170 


5 


.36 


1 


.58 





,64 


1 


.35 


-1.78 


-0 


,52 


— 1 


.61 


—2 


.80 


-2 


.84 


-2.84 


2 


NGC3599 


21.7 


80 


13 


.01 


1 


.66 





.79 





.50 


— 1.64 


— 1 


,01 


— 1 


.65 


— 1 


.34 


— 1 


.35 


-1.34 


1 


NGC3605 


21.7 


103 


9 


.14 


1 


.26 





67 





.92 


— 1.82 


-0 


75 


— 1 


.64 


-1 


.84 


-1 


.86 


-1.85 


1 


NGC3608 


21.7 


195 


1 


.05 


1 


,33 





,00 


1 


.88 


— 2.32 


-0 


29 


— 1 


.53 


—2 


.80 


—2 


.90 


— 2.91 


1 




13.6 


195 





.72 


1 


.58 





00 


1 


.88 


— 2.28 


-0 


.37 


-1 


.63 


-3 


.03 


—3 


.15 


-3.16 


3V 




13.6 


195 





.78 


1 


.57 





,00 


1 


.92 


— 2.32 


-0 


39 


— 1 


.67 


—3 


.11 


-3 


.25 


-3.26 


31 


NGC3900 


29.4 


118 





.29 


1 


.66 





,51 


1 


.39 


-2.24 


-0 


,75 


—2 


.37 


-3 


.65 


-6 


.00 


-6.00 


2 


NGC4026 


17.0 


195 





.88 


1 


50 





,68 


1 


.51 


-1.92 


-0 


,41 


-1 


.84 


-3 


.41 


-3 


.94 


-4.07 


2 


NGC4143 


17.0 


270 


1 


.26 


2 


.18 





.59 


1 


.20 


-1.35 


-0 


,26 


-1 


.45 


-3 


.23 


—3 


.35 


—3.42 


2 


NGC4150 


9.7 


85 


1 


.23 


1 


.67 





,58 


1 


.18 


-2.28 


-1 


,09 


—2 


.20 


-2 


.69 


—2 


.79 


-2.74 


2 


NGC4168 


38.8 


185 





.95 


1 


,50 





14 





.76 


— 1.11 


-0 


10 


— 


.82 


— 1 


.24 


-1 


.25 


— 1.25 


1 


NGC4239 


16.3 


60 


14 


.53 





.96 





65 





.39 


-1.78 


— 1 


,09 


— 1 


.71 


— 


.75 


— 


.76 


-0.75 


1 


NGC4261 


35.1 


326 


2 


.38 


1 


43 





00 


1 


.43 


— 1.32 





.15 


— 


.86 


-2 


.38 


—2 


.41 


— 2.44 


2 


NGC4278 


9.7 


250 


1 


.63 


1 


.39 





02 


1 


,96 


— 2.10 


— 


16 


— 1 


.50 


-3 


.25 


-3 


.36 


-3.40 


2 




9.7 


250 


1 


45 


1 


,32 





00 


1 


,98 


— 2.13 


— 


09 


— 1 


43 


—3 


.12 


-3 


.23 


-3.27 


3V 




9.7 


250 


1 


.25 


1 


.46 





00 


1 


,90 


— 2.04 


— 


15 


— 1 


.43 


-3 


.13 


-3 


.22 


-3.25 


31 


NGC4291 


29.4 


278 


2 


07 


1 


.48 





,02 


1 


,89 


-1.92 


-0 


.14 


— 1 


.44 


-3 


.28 


-3 


.38 


-3.43 


2 


NGC4365 


23.5 


262 


2 


.06 


1 


.27 





15 


1 


,43 


-1.51 





07 


— 


.96 


-2 


.25 


-2 


.28 


-2.30 


1 




13.8 


262 


1 


.67 


1 


,46 





,11 


1 


.58 


-1.63 


-0 


,04 


— 1 


.15 


-2 


.66 


-2 


.71 


-2.73 


3V 




13.8 


262 


1 


.52 


1 


.49 





04 


1 


.61 


-1.67 


-0 


06 


— 1 


.17 


-2 


.69 


-2 


.73 


-2.75 


31 


NGC4374 


16.8 


296 


2 


.15 


1 


.50 





.13 


1 


.53 


-1.47 





03 


-1 


.05 


-2 


.66 


—2 


.70 


-2.72 


2 


NGC4387 


16.3 


105 


3 


.36 


1 


.59 





72 





.38 


— 1.24 


-0 


66 


-1 


.33 


— 1 


.24 


-1 


.26 


— 1.25 


1 


NGC4406 


16.8 


250 


3 


.31 


1 


.16 





.00 


1 


.70 


— 1.84 





.03 


— 1 


.19 


—2 


.61 


-2 


.67 


-2.69 


2 




13.8 


250 


4 


.13 


1 


.05 





04 


1 


,80 


-1.95 





.13 


— 1 


.12 


—2 


.51 


—2 


.58 


— 2.61 


3V 




13.8 


250 


3 


.32 


1 


.07 





00 


1 


,81 


-1.97 





12 


-1 


.14 


-2 


.55 


—2 


.63 


-2.66 


31 


NGC4417 


16.8 


84 





.87 


1 


.77 





,71 





,10 


-1.28 


-0 


,90 


-1 


.80 


-1 


.75 


— 1 


.80 


-1.78 


2 


NGC4434 


16.3 


115 





.98 


1 


.78 





,70 





,60 


-1.51 


-0 


.74 


-1 


.77 


-2 


.31 


-2 


.39 


-2.35 


1 


NGC4458 


16.3 


105 


5 


.26 


1 


.43 





49 


2 


.02 


— 2.81 


— 


.99 


— 2 


.51 


—3 


.87 


—4 


.04 


—3.96 


1 


NGC4464 


16.3 


125 


1 


.64 


1 


.68 





,88 





.96 


-1.81 


-0 


,81 


-2 


.08 


-3 


.11 


-6 


,00 


-3.35 


1 


NGC4467 


16.3 


72 


7 


.52 


2 


.13 





,98 





.02 


-1.38 


-1 


.10 


-2 


.75 


-3 


.38 


-6 


.00 


-6.00 


1 


NGC4472 


16.3 


300 


2 


.08 


1 


,17 





04 


1 


.51 


-1.51 





,24 


-0 


.82 


-2 


.21 


-2 


.24 


-2.26 


1 




16.8 


303 


1 


.89 


1 


,29 





04 


1 


.48 


-1.42 





,20 


-0 


.84 


-2 


,29 


-2 


.32 


-2.34 


2 


NGC4478 


16.3 


135 


3 


.32 





.84 





,43 


1 


.75 


-2.66 


-0 


24 


-1 


.17 


-1 


.39 


-1 


,42 


-1.42 


1 


NGC4486 


16.3 


360 


2 


.82 


1 


.39 





,25 


1 


.32 


-1.01 





,36 


-0 


.63 


-2 


.28 


-2 


.29 


-2.31 


1 


NGC4486B 


16.3 


200 


2 


.78 


1 


.33 





14 


2 


.24 


-2.63 


-0 


,50 


-2 


.10 


-3 


.90 


-4 


.33 


-4.38 


1 


NGC4551 


16.3 


121 


2 


.94 


1 


.23 





,80 





.54 


-1.30 


-0 


,52 


-1 


.36 


-1 


.56 


-1 


.59 


-1.58 


1 


NGC4552 


16.3 


260 


1 


.48 


1 


.30 





00 


1 


.93 


-2.08 


-0 


,05 


-1 


.36 


-3 


.02 


-3 


.12 


-3.16 


1 




13.8 


260 


2 


.17 


1 


.06 





,00 


2 


.06 


-2.25 





,14 


-1 


.23 


-2 


.79 


-2 


.92 


-2.98 


3V 




13.8 


260 


2 


.10 


1 


.08 





04 


2 


,06 


-2.23 





,12 


-1 


.27 


-2 


.87 


-3 


.00 


-3.06 


31 


NGC4564 


16.3 


165 





25 


1 


.90 





,05 


1 


.61 


-2.17 


-0 


,45 


-1 


.80 


-3 


.02 


-3 


.33 


-3.29 


1 


NGC4570 


16.3 


195 


3 


.72 


1 


.49 





,85 


1 


.04 


-1.44 


-0 


,36 


-1 


.43 


-2 


.57 


-2 


.65 


-2.65 


1 


NGC4589 


30.0 


228 


1 


.09 


1 


.18 





,11 


1 


.93 


-2.27 


-0 


,07 


-1 


.32 


-2 


.64 


-2 


.75 


-2.77 


2 




24.4 


228 





.43 


1 


.62 





,00 


1 


.55 


-1.80 


-0 


05 


-1 


.12 


-2 


.30 


-2 


.35 


-2.36 


3V 




24.4 


228 





.50 


1 


.58 





,01 


1 


.68 


-1.95 


-0 


11 


-1 


.24 


-2 


,50 


-2 


.58 


-2.59 


31 


NGC4621 


16.3 


250 





.19 


1 


.71 





,50 


1 


.16 


-1.42 


-0 


.10 


-1 


.90 


-3 


.94 


-6 


.00 


-6.00 


1 


NGC4636 


16.3 


210 


1 


.64 


1 


,33 





,13 


1 


.14 


-1.38 


-0 


,04 


-0 


.90 


-1 


.77 


-1 


.78 


-1.78 


1 




17.0 


207 


1 


.69 


1 


.56 





13 


1 


.07 


-1.30 


-0 


14 


-0 


.96 


-1 


.84 


-1 


.85 


-1.85 


2 


NGC4649 


16.3 


360 


2 


.00 


1 


.30 





,15 


1 


.55 


-1.34 





,32 


-0 


.79 


-2 


.50 


-2 


.54 


-2.57 


1 


NGC4697 


11.2 


175 


24 


.86 


1 


.04 





,74 


1 


.27 


-1.64 


-0 


,11 


-1 


.26 


-2 


.24 


-2 


.28 


-2.28 


1 


NGC4742 


13.3 


105 


48 


.60 


1 


.99 


1 


,09 





.78 


-1.83 


-1 


,05 


-6 


.00 


-6 


.00 


-6 


.00 


-6.00 


1 
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Table 1 — Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Name 


D 


a 


a 


P 


7 


log K b 


log r b 


log -Rein 


log -R ra d 


log(Mcoro) 


iog(Mcoro) 


log(Mcorc) 


Ref 




(Mpc) 


(km/s) 










(arcsec) 


(arcsec) 


(arcsec) 










NGC4874 


99.5 


290 


2.33 


1.37 


0.13 


0.71 


-0.68 


0.25 


-0.39 


-1.08 


-1.08 


-1.09 


1 


NGC4889 


99.5 


350 


2.61 


1.35 


0.05 


1.06 


-0.88 


0.33 


-0.47 


-1.66 


-1.68 


-1.69 


1 


NGC5273 


21.3 


52 


7.03 


1.32 


0.37 


0.48 


-1.93 


-1.26 


-1.79 


-0.72 


-0.73 


-0.73 


2 


NGC5813 


30.2 


225 


2.15 


1.33 


0.08 


1.50 


-1.72 


-0.13 


-1.19 


-2.42 


-2.46 


-2.47 


1 




21.2 


225 


1.77 


1.41 


0.03 


1.57 


-1.79 


-0.17 


-1.27 


-2.58 


-2.62 


-2.63 


3V 




21.2 


225 


1.67 


1.46 


0.01 


1.58 


-1.80 


-0.20 


-1.29 


-2.62 


-2.67 


-2.68 


31 


NGC5838 


28.5 


290 


2.57 


1.87 


0.93 


0.83 


-0.98 


-0.12 


-1.29 


-2.97 


-6.00 


-6.00 


2 


NGC5845 


30.1 


260 


1.27 


2.74 


0.51 


1.03 


-1.27 


-0.35 


-1.44 


-3.14 


-3.23 


-3.27 


1 


NGC5982 


38.7 


256 


1.73 


1.28 


0.06 


1.65 


-1.80 


-0.02 


-1.16 


-2.57 


-2.62 


-2.65 


2 




39.9 


256 


2.15 


1.19 


0.12 


1.64 


-1.79 


0.03 


-1.13 


-2.50 


-2.55 


-2.58 


3V 




39.9 


256 


2.17 


1.19 


0.12 


1.64 


-1.80 


0.03 


-1.13 


-2.51 


-2.57 


-2.59 


31 


NGC6166 


120.0 


300 


3.32 


0.99 


0.08 


0.67 


-0.68 


0.35 


-0.27 


-0.74 


-0.75 


-0.76 


1 


NGC6340 


22.0 


137 


2.46 


1.28 


0.59 


1.61 


-2.28 


-0.62 


-1.91 


-2.97 


-3.09 


-3.06 


2 




18.0 


137 


1.73 


1.24 


0.72 


1.38 


-2.08 


-0.57 


-1.80 


-2.70 


-2.83 


-2.79 


4 


NGC7332 


21.7 


130 


4.25 


1.34 


0.90 


1.12 


-1.88 


-0.68 


-1.87 


-2.76 


-2.93 


-2.86 


1 


NGC7457 


12.3 


77 


2.32 


1.03 


0.35 


1.22 


-2.46 


-0.90 


-1.83 


-1.67 


-1.71 


-1.69 


2 


NGC7626 


45.3 


273 


1.84 


1.30 


0.36 


1.57 


-1.64 


0.01 


-1.15 


-2.65 


-2.71 


-2.75 


2 




49.5 


273 


1.53 


1.23 


0.00 


1.64 


-1.77 


0.08 


-1.03 


-2.40 


-2.45 


-2.47 


3V 




49.5 


273 


1.29 


1.27 


0.00 


1.65 


-1.78 


0.07 


-1.04 


-2.42 


-2.47 


-2.49 


31 


NGC7743 


24.4 


85 


5.36 


1.38 


0.50 


1.47 


-2.48 


-1.04 


-2.19 


-2.76 


-2.82 


-2.79 


2 


NGC7768 


110.0 


290 


1.92 


1.21 


0.00 


1.39 


-1.46 


0.17 


-0.82 


-2.04 


-2.07 


-2.08 


1 



Note. — Col. 2: Distance (assuming Hq — 75 km s _1 Mpc -1 ). Col. 3: Central stellar velocity dispersion. Cols. 4-6: Fitted Nuker law parameters. 
Cols. 7-10: Logarithms of the lensing strength Kb, break radius r^. Einstein radius i? c i n , and radial critical radius i? ra d ■ Note that if any quantity is zero, 
we arbitrarily set its logarithm to —6. Col. 11: Logarithm of mean core image magnification for the Nuker law lens. Cols. 12-13: Logarithm of mean core 
image magnification for the Nuker law lens plus a supcrmassivc black hole normalized by the M,-<r correlations of Gebhardt et al. (2000) and Merritt & 
Fcrrarese (2001), respectively. Col. 14: References as follows: (1) Faber ct al. (1997): (2) Ravindranath et al. (2001); (3V) and (31) V and I band samples, 
respectively, from Carollo ct al. (1997); (4) Carollo & Stiavelli (1998). 



